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■ Abstract 

We consider the A'^^-flavour Schwinger Model on a thermal cylinder of circumference 
(3 = 1/T and of finite spatial length L. On the boundaries =0 and = L the fields 
are subject to an element of a one-dimensional class of bag-inspired boundary conditions 
which depend on a real parameter 9 and break the axial flavour symmetry. For the cases 
Nf = 1 and Nf — 2 all integrals can be performed analytically. While general theorems 
do not allow for a nonzero critical temperature, the model is found to exhibit a quasi- 
phase-structure: For finite L the condensate - seen as a function of log(T) - stays almost 
constant up to a certain temperature (which depends on L), where it shows a sharp 
crossover to a value which is exponentially close to zero. In the limit L —t oo the known 
behaviour for the one-flavour Schwinger model is reproduced. In case of two flavours 
direct pictorial evidence is given that the theory undergoes a phase-transition at Tc = 0. 
The latter is confirmed - as predicted by Smilga and Verbaarschot - to be of second order 
but for the critical exponent 5 the numerical value is found to be 2 which is at variance 
with their bosonization-rule based prediction (5 = 3. 
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Figure 1: The (dimensionless) condensate as a function of log{kT / fii) at fixed 

box-length L = 1/ fii (where /ii = |e|/7r^/^ is the induced single-flavour Schwingermass (^)) 
for Nf = 1 and Nf = 2. 



1 Introduction 

For realistic gauge field theories like QCD(4) it is in general an unsolved problem to determine 
their phase-structure (e.g. as a function of the fermion masses rrii . . .m^j:) analytically. For 
this reason, one may either determine their phase-structure approximatively or try to attack 
the question in some simpler models analytically |2| . 

For many questions arizing in QCD, the Schwinger model (QED in two dimensions 
with one or A'^^ massless fermions) has proven to be an interesting testing ground: It has shed 
some light on such longstanding problems of QCD as the C/(l)yi-problem |^ and it has been 
used to test the validity of the Instanton Liquid Picture |5|, ^ . 

However, for investigating symmetry-breakdown any two-dimensional model field theory 
doesn't seem to be of interest: There is no spontaneous breaking of a continuous global 
symmetry with associated Goldstone bosons in two dimensions Q and a symmetry which is 
anomalously broken can't get restored at finite temperature Q. 

It is the aim of the present paper to show that - in spite of the truth of this conventional 
wisdom - the Schwinger model exhibits an interesting quasi-phase-structure: The chiral con- 
densate ib75)'i/') (which is used to probe the chiral symmetry) shows - as a function of 
the log of the temperature - a sharp crossover-behaviour: For any finite box-length L there is 
a well defined low-temperature regime where the condensate stays almost constant (the value 
depends on L) and there is a "critical temperature" where the condensate decays (through 
a fairly well localized symmetry-quasi-restoration process) to a value which is exponentially 
close (but not equal) to zero. 

Moreover, as the length L is sent to infinity, this behaviour is shown to provide direct 
evidence that the two-flavour Schwinger model exhibits a phase-transition at Tc = and that 
the latter is of second order - which is an assumption-free rederivation of a recent claim by 
Smilga and Verbaarschot |^. 

The Schwinger model has had a great impact on the development of field-theoretic ideas 
and techniques: The original quantization on the plane suffered from the deficit that 
a direct calculation of the condensates (V'^V') ^iid gave vanishing results whereas 



an indirect determination via the clustering theorem led to the standard nonzero value |11|. 
One decade ago the Schwinger model has been quantized on compact manifolds without 
boundaries, the first ones being the (euclidean) sphere [f^] and the (euclidean) torus 
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|13|. In either case direct calculations for the chiral condensates (ip^ip) ^-iid ('0^75 V') were 
found to yield nonvanishing results for a finite volume of the manifold. 

Here we shall consider the model on a finite-temperature-cylinder with (anti-)periodic 
boundary conditions in the euclidean timelike direction, but with noncyclic boundary condi- 
tions in the spatial direction: On the two spatial ends (at x^ = and x^ = L) some chirality- 
breaking (XB-) boundary-conditions are imposed. 

The motivation to use these XB-boundary-conditions stems again from QCD. There one 
would like to give a proof (without involving any assumption) that in the chiral limit the axial 
flavour symmetry SU^Nf)^ is spontaneously broken. The natural way to study a symmetry 
which is expected to be broken spontaneously is to break it explicitly and to try to determine 
how the system behaves in the limit when the external trigger is softly removed. For both 
QCD and the A^-flavour Schwinger model this means that one has to break the axial A^- 
flavour symmetry and then try to determine how some observables which are sensitive to its 
breaking do behave in the limit where the trigger term is removed. There is a long series 
of efforts in the literature which try to achieve this goal by studying either QCD or the A^- 
flavour Schwinger model with a small fermion mass. In this approach the task is to determine 
how the chiral correlators behave in the limit where the fermion masses tend to be tiny as 
compared to the intrinsic mass of the theory. The problem with this particular method of 
breaking the axial flavour symmetry is that the value of the chiral condensate is related to 
the mean level density of the eigenvalues of the Dirac operator in the infrared |14], but the 
spectral density of the massive Dirac operator in a given gauge-field background is not (yet) 
known in general. 

Together with A.Wipf, we have previously explored the alternative of breaking the chi- 
ral symmetry through introducing boundary conditions for the fermions rather than giving 
them a mass [15|. There we dealt with euclidean U{Nc) and SU{Nc) gauge-theories with 

massless flavours quantized inside an even(d = 2n)-dimensional ball Bj^ with boundary 
on which the XB-boundary-conditions considered in [^] have been imposed. These 
boundary conditions relate the different spin-components of each flavour on the boundary 
and are neutral with respect to vector-fiavour-transformations. As a consequence the (gauge- 
invariant) fermionic determinant is the same for all fiavours. The most important result of 
this work was the observation that the XB-boundary-conditions proved to be equally suited 
to trigger a chiral condensate as small fermion masses are. From a technical point of view, the 
XB-boundary-conditions turned out to be very convenient: In abelian theories, the (mass- 
less) fermionic determinant can be calculated in an arbitrary gauge-field background (the 
fields subject to the XB-boundary-conditions) and thus we derived in two dimensions the 
analytical expressions for the chiral correlators (exploiting the rotational symmetry of B"^). 

In a second work we examined whether the approach of breaking the SU^Nf)^- 
symmetry by boundary-conditions can be extended to gauge-systems at finite temperature. 
Choosing the Schwinger model as a testbed the answer was in the affirmative: First of all we 
found that most of the nice features associated to the quantization on contractible manifolds 
(topologically equivalent to B'^^) persist in the case of the finite temperature cylinder with 
XB-boundary-conditions at the spatial ends: The configuration space is still topologically 
trivial (i.e. without disconnected instanton sectors) and in particular there are no fermionic 
zero-modes (which usually tend to complicate the quantization considerably |jl3|). The tech- 
nical difficulty we had to deal with is the fact that on non-contractible manifolds (e.g. a 
cylinder) the standard decomposition-technique for the Dirac operator, on which the quanti- 



zation adopted in [15[| heavily relied, could no longer be used. For this reason |17[ turned out 
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to be a rather technical paper, mostly devoted to show how this difficulty can be overcome. 

The present paper concentrates on the physics which can be addressed from this setup. 
To be definite: The Schwinger model 

S[A,i:\^] = Sb[A] + SF[A,i:\^] 

Nf (1) 

Sb = IJf^uF^u , Sf = I i'iipi^n 

M n=lM 

in d = 2 dimensions is studied on the manifold 

M = [0, /3] X [0, L] 3 (2) 

with volume V = PL. In euclidean time direction, the fields A and ip are periodic and 
antiperiodic respectively with period p. Hence = and = (3 are identified (up to an 
eventual minus sign) and the manifold is a cylinder. At the spatial ends of the cylinder (at 

= and x^ = L) specific XB-boundary-conditions (to be discussed below) are imposed. 
Through explicit numerical evaluation of the analytical formulae we demonstrate that the 
model exhibits a quasi-phase-structure as described above. In particular we find that the 
two-flavour Schwinger model (limit L oo) exhibits a (true) phase-transition at Tc = and 
that the latter is of second order - which is a rederivation of a result which was obtained by 
Smilga and Verbaarschot |^ using completely different technical tools. 

For notational simplicity we use the abbreviations 

"=2l ' ^ = T 
as well as the dimensionless inverse temperature and box-length 



which are built from /3, L by rescaling them with the Compton-wavelength associated to the 
single-flavour Schwinger-mass /xi where 



(5) 



is the generalization of /ii to the case of A^-flavours. 

This paper is organized as follows: In section 2 we briefly review the quantization of 
the model subject to the XB-boundary-conditions. In section 3 the analytical results for 
the chiral condensate (^/'^^^(l ib75)^/')(x) are re-expressed using theta- functions. In section 
4 the remaining (c-number-)integrals within the analytic formulae for arbitrary Nf are fully 
performed for the cases Nj = 1 and Nf = 2. This allows us to show explicit numerical 
evaluations of our analytical findings for these two cases revealing the two quasi-phases as 
described above and resulting in an explicit plot of the "quasi-phase-structure" figure 0. In 
the concluding section 5 we try to illustrate our results and to relate them to the work by 
Smilga and Verbaarschot |0] and others. 
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2 Quantization with Thermal and XB-Boundary-Conditions 



Here we shall give a short review of the conceptually relevant aspects of the quantization. For 
technical aspects the reader is referred to [p^] . 

2.1 Chirality Breaking Boundary Conditions 

The chirality-breaking (XB-)boundary-conditions as discussed in [16| can be motivated by 
the request that the Dirac operator iJp is symmetric under the scalar product (x, V') '■= 
J x^i^ dPx, which leads to the condition that the surface integral i § x^lni^ ds vanishes, where 
In = (7i = n^^fj_ = 1^ and is the outward oriented normal vectorfield on the boundary. 
Imposing local linear boundary conditions which ensure this requirement amounts to have 
X^lni^ = on the boundary for each pair. A sufficient condition is to have all modes obeying 
ij) = Bip on the boundary, where the boundary operator B (which is understood to act as the 
identity in flavour space) is required to satisfy B'^'jnB = — 7n and B^ = 1. We shall choose 
the one-parameter family of boundary operators |16| 

B = Be:= i^^e'^'^ln (6) 

which break the 75-invariance of the theory, thus making the Nf flavour theory invariant under 
SU{Nf)y instead of SU{Nf)L x SU^Nf)^. They will be supplemented by suitable boundary 
conditions for the gauge-field. These boundary conditions prevent the C/(l)-current from 
leaking through the boundary as they ensure j-n = ip^jri'ip = on the boundary. 

For explicit calculations we shall choose the chiral representation 70 = <ti,7i = (T2 and 
75 = (J3 which implies that the explicit expressions for the boundary operators at the two 
ends of the cylinder (0) take the simple form 



Bl = -{ -en (atxi=0) and Br = + { (atxi=L) . (7) 





2.2 Immediate Consequences for the Spectrum 

The decision to quantize with boundary condition ip = Bgip with the boundary operator (jBI) 



has immediate consequences [15, 17, 18 



(i) The Dirac operator has a discrete real spectrum which is asymmetric w.r.t. zero. 

(ii) The spectrum is empty at zero, i.e. the Dirac operator has no zero modes. 

(Hi) The instanton number q = ^ e^iuF^iu = ^ / ^ ^ is not quantized. 

The first property already indicates that we are not in the situation covered by the Atiyah- 
Patodi-Singer-index-theorem. The second property implies that the generating functional for 
the fermions in a given gauge-field background A 

ZF[A,rilri] = f D-^^D^ e-I^^'^+'I^^''-'I''^'f' (8) 
Np J 

is indeed given by the textbook formula 

Z^\A v] - '^^^'^'^^ eI^'^'V>)~'^ (9) 
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and the chiral expectation values follow by taking the logarithmic derivative 



Throughout is the free parameter in the boundary operator (j^. The fact that the 
Feynman-Hellmann-boundary-formula ^Afc = —\k{'^kilb'^k) HI] (where the denote the 
eigenvalues of ilp) still holds true on the cylinder was the basis for the analytic determination 



of the ^-dependence of the fermionic determinant det£i(i^) performed in |17] 



2.3 Neither Integer nor Fractional but Real Instanton Number 

On a cylinder of finite spatial length the decomposition of a gauge potential is |17| 

eAi= +do^ + dix 

where (j) obeys Dirichlet boundary conditions at the ends (x^ = 0, L) and x is a pure gauge 
degree of freedom which fulfills x(0) = and c G [—1/2, l/2[ is the constant harmonic 

part. Thus the Dirac operator il/)= i^p,{d^ — ieA^) may be factorized 

ilp= GUpoG (12) 

where ipQ = ^^{id^ + 2ttc/ j3) + 7^i3i is the Dirac operator with the scalar parts switched 
off and G = (Ma.g{g*~^ , g) contains the prepotential g := e~^'^^'^^^ which is an element of the 
complexified gauge-group U{1)* = x R+. 

On the cylinder there is a one-to-one-correspondence between (p and e-Foi if obeys 
Dirichlet boundary-conditions at the two ends. The general field (j) may be decomposed as 

,4- ,27rmj;°, . ,7rnx\ ,_ . Inmx^ . . irnx^ 
(P = 2^ 2^ <l>rnn cos( — - — ) sm(— ^) + <f>^^ sm( — - — ) sm(— ^) (13) 

m>On>l '^^ ' 

with coefficients (pmn ^ ^ decaying rapidly enough to make the series converge. The instanton- 
number is given by 

n>l 



where the decomposition (11) and the expansion ( |l3|) have been used. Given the result ( [14| 



it is obvious that the instanton-number q may take any real number. 
2.4 Fermionic Propagator w.r.t. Boundary Conditions 

In order to calculate the condensates one needs the Green's function Sq of the Dirac operator 
Upon the cylinder subject to the XB-boundary-conditions. In addition to the defining relation 
{ip Sg){x,y) = 5{x — y) this Green's function obeys the boundary-conditions 

Se{x''+P,x\y) = -Se{x,y) (15) 
{BLSe){x\x^ = Q,y) = Seix"" ,x^ = Q,y) (16) 
{BRSe){x\x^=L,y) = Seix"* ,x^ = L,y) (17) 
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with -B^/ij defined in (0) plus the adjoint relations with respect to y. The dependence of the 
gauge-potential has not been made explicit, since from the factorization-property ([T^ ) for the 
Dirac-operator it follows at once that Sg is related to the Green's function Sg of ipo as 

Se{x,y) = G-\x)Se{x,y)G^-\y) . (18) 

Since the field (j) obeys Dirichlet boundary-conditions at the ends of the cylinder, the de- 
formation matrix G appearing in (12) is unitary there and the boundary-conditions ([T5|-p!^) 
transform into the identical ones for 

Se{x,y) -- 

where the indices refer to chirality. This Green's function - which carries the full c-dependence 



of Sg{x,y) ~ has been determined analytically [|17[ to read 



Sg{x, y) = — - (-1)^"^+"^ • e2--((-°-^°)//^-") . ( _f/!>" [y/_"™^ ) , (19) 

where r„m = {x^—y^) + i{x^+y^) — {nP + 2imL) and Snm = {x^—y^) + i{x^—y^) — {nP + 2imL). 

From ([l9| ) it follows that the +-|- and elements at coinciding points inside the cylinder 

take the forms 

sg{x,x)^ = ±^E(-r ■ . 

_l_e^^ cos(27rnc) sin(7r^)ch(7rnT) — sin(27rnc) cos(7r^)sh(7rnr) 

1l sin2(7re) +sh2(7rnr) ^ ' 

p±e -27r(m+Oc/T 

~Se[x,x)^ = ±_ ^ (-)™— 

n^z smh(^(m + i)/T) 

I V r + ^^"^/^^ ~ sh(27r(m + g)c/r) 

both valid for c G [— ^, The two forms (^, |2l|) are equivalent but enjoy good convergence 
properties in the two regimes j3 ^ L and P <^ L respectively. From (|l^) one sees that the 
chirality violating entries of the fermionic Green's function lie on the diagonal and take the 
form 

Se{x; x)^ = e^20(-)5^(2;; x)^ (22) 



where S±± plugged in from (2C, |2l] ) depends only on the harmonic part c in the decomposition 
([ll|) of the gauge-potential. 



2.5 Fermionic Determinant w.r.t. Boundary Conditions 



The arduous step is the computation of the ^-dependence of the fermionic-determinant |17|. 
The Dirac-operator and the boundary-conditions are both flavour-neutral. Thus the deter- 
minant is the same for all flavours and it is sufficient to calculate it for one flavour. For 
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the explicit calculations we used the gauge-invariant (^-function definition of the determinant 

logdete(i^):=^logdete(-^):=-^^ _Ce{-f,s) (23) 

and calculated the 0-dependence of the ^-function by means of a boundary-Feynman-Hell- 
mann-formula. Denoting {^k\k G N} the (positive) eigenvalues of — the corresponding 
("-function is defined and rewritten as a Mellin-transform in the usual way 

oo 

C,(.) G(-#,.) = fZ^/ t^-nr,(e-*(-^)) dt (24) 

for Re(s) > d/2 = 1 and its analytic continuation to Re(s) < 1. 
The general task was to compute the normalized determinant 



dete(i^) _ dete{iiPi,c 



deto(i^) deto(i-^o,o) 



(25) 



where the first suffix on the r.h.s. indicates whether the scalar part is switched on and the c 
refers to the harmonic part in the decomposition (0). Together with 6 we thus have three 
parameters to switch off and this leaves us with 3! = 6 possible choices how to compute the 
functional determinant (^5|) in terms of three factors where each involves one switching only. 
We explicitly followed two of the six choices and found them to agree. 

Our first choice was to calculate the functional determinant according to 



deie{ilpi,c) _ dete(z^i,J deto(i^i,c) deto(i^ 



deto(i^o,o) deto(i^i,c) deto(i-^o,c) deto(i^o,( 
where we got for the first factor the explicit expression 

dete(i^i,c) 



(26) 



deto(i^i,c! 
to be multiplied with the result 

deto(z-^i,e) 
deto(i-^o,c) 



exp{-— j ee^,F^,} = exp{-— J A0}. (27) 



for the second factor. It is worth noting that with the first choice the term linear in 9 in the 
effective action stems from a volume-term (i.e. oi(-)) in the Seeley-DeWitt-expansion. 
Our second choice was to calculate the functional determinant according to 

dete(i^i,c) ^ dete(z^i^c) dete(i^o,c) deto(i^o,c) .^g) 
deto(i^o,o) det0(i^o,c) deto(i-^o,c) deto(i^o,o) 

where we found for the first factor the explicit expression 

dete(i^o,c) ^^2tt J ^ ^ 27r J ^ ^ 

to be multiplied with the result 

dete(«^o,c 



exp{0} = 1 (31) 
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for the second factor. It is worth noting that with the second choice the term hnear in 9 in 
the effective action stems from a boundary-term (i.e. 6i(.)) in the Seeley-DeWitt-expansion. 

In summary, the scattering of the fermions off the boundary generates a CP-odd term 
linear in 9 in the effective action for the gauge-bosons which may be seen as a two-dimensional 
artificial analogue of the QCD-0-term. 

The remaining task (which is to calculate the common third factor in the factorizations 
(p6|) and (2£) of the functional determinant) was addressed by rewriting its logarithm as 



deto(z-^o,cj I [ d d ^ 
deto(«-^o,o) 2 J ds s=o"c 

and constructing the s-derivative at s = of the c-derivative of Co(— -^oc explicitly. For 
that aim we computed the heat-kernel of the operator 



mc = -{ido-27ric/(3r + dtl2 



for 9 = on the half-cylinder (see appendix of [|17|) and from this expression we were able to 
derive 

r(c) . -log ^!^4f^ = ii:Vi)"-^" "y!7' m 

deto(«-^o,o) vr ^ [nPy + (2mLY 



as the result for the negative of the logarithm of the last factor of ( |26D and (^) . In (|3^) and in 
the following the prime in the sum denotes the omission of the contribution from m = n = 0. 
By performing either the sum over m or the sum over n in (|33|) one gets 



r(c) = y tH. ^o^'^^^c) - 1 (34) 
^' 4l n sh(n^/?/2L) ^'^^ 

- (-^r ch(4m7rcL//3) - 1 2nL , 

- M2m.L/(3) + (3 ' ^^^^ 

both valid for c G [—1/2,1/2] and periodically continued otherwise. These two equivalent 
forms will be useful in the low- and high- temperature expansion of the condensates. 

2.6 Effective Action 



The final step is to combine the classical (euclidean) action of the photon field, rewritten in 
the variables ( [ll| ) 

with the result for the functional determinant (p5|). Collecting the contributions (27, |28| ) or 



Sb[<P] = -F^,F^, = -^A^A<I) (36) 



(|30|, |3l|) as well as (|34D or (|35|) and adding the classical action (^) one ends up with the 



effective action (which, of course, does not contain the gauge degree of freedom x) 

r = Te,N, [c, ct>]=Nf r(c) + Te,N, [4>] (37) 
where r(c) has been given in ( p4| , |3^ ) and Tg^j^^lcp] is 

^e,Nf M = ^ { / '^^''^ - (t>^<P + ■ 1 A</.| (38) 
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where the fact that the functional determinant is the same for all flavours has been used. In 
( |38| ) /i is the Schwinger mass (||) which is the analog of the j^'-mass in 3-flavour-QCD. 
In summary, the functional measure takes the form 

df,e[A] = ^ e-'^^'^'f^''^^ dc D4> <5(x) Dx (39) 

where we have taken into account that the gauge-variation of the Lorentz gauge-condition 
F := dfj^A^ = Ax and the Jacobian of the transformation from {A} to the variables {(p, c, x} 
are independent of the fields. Actually, the corresponding determinants cancel each other. 

We conclude that the expectation- value of any gauge-invariant operator O (which will not 
depend on x) is given by 

(O) = /^^^^ Oe-'-^A-^^ 



with Tg^Nf [c, 4>] given by (0, and (|3|, ^ . 

3 Condensates in a Finite Volume 

The general result (^) may be applied to calculate the chiral condensates 

{ip^{x)P±ij{x)) = _r, „ [c ^1 (41) 



with 5*61 from (22) and Tg from (37). Both the (exponentiated) action and the Green's function 
factorize into parts which only depend on c and (/>, respectively. Thus (^) factorizes as 

{ipHx)P±ij{x)) = C^{x) ■ D^{x) (42) 

with x'^-independent factors 

which depend on the parameters 6,Nf,P,L. Here and below the c-integrals extend over the 
period [—1/2, 1/2], whereas the field (p is subject to Dirichlet boundary conditions at the two 
ends = and x^ = L. The next step is to evaluate the factors and in ( ^2|) for given 
circumference /3 and length L of the cylinder and given values of 9 and Nf. 

3.1 Harmonic Integral 

In order to evaluate the first factor in ( p2|) it is worth noticing that the two forms (^, ^5|) 
allow one to write the factor exp{— r(c)} in the two equivalent versions 

r(c) ^ OsicJr) 
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^3(0, i/r) ^ ^ 

Here we employed the notation 

03(u, uj) = Y, e2™"g"' = 1 + 2 ^ cos(2n7ru)g"' (g = e*^"^) (47) 

neZ n>l 

for the parameters oo = it and u) = i/r (giving real nome q g]0, 1[) respectively. In deriving 
( |45| ) and (|4^) we used the infinite-product expansion [^] 

^3(^,0;)= + V""^cos(27rn) + ^2(2n-i)) (4g) 

n>l 

and the addition theorem [E3I 



In ( ^^i^+^'^i ) = 4 y izil::^— sin(2^nn) sin(2vrnt;) (49) 
respectively. The Poisson resummation lemma makes sure that 



^ g-(x+n)2.t ^ ^^/^ ^ n Ag2™x (x e iZ, t > 0) (50) 

from which we derive the theta-function duality relation 

^e^{c,iT)=e-^''^l^ e^{±ic/T,i/T) {ceR,T>0) (51) 

which in turn allows to confirm the identity of the right-hand sides of (H^) and (46). 



Starting from expression (|43| ) and plugging in the Green's function ( poj) as well as (^) or 
alternatively the Green's function (^) along with (^) one ends up with 

^ - \ ^^^ ^\n sin«)ch(7rnr) / cos(27rnc) 6'f^'(c, ir) dc 

^ " ^ 4L sin2(vrO +sh2(vrnr) ' / Q^,^ {c,^r) dc ^ ^ 



/ch(27r(m+^)c/r)e-^/^^ e3^ (ic/r, i/r) dc 



m^z sh(7r(m+e)/r) j g-7V^^cVr 6/f^(ic/r,i/r) dc 

(53) 

which is independent of as required by translation-invariance. In the last step we have taken 
advantage from the fact that ^3(., .) is symmetric in its first argument and the c-integration 
is from -1/2 to 1/2. 

3.2 Scalar Integral 



For the evaluation of the second factor in (42) we recall that the integration extends over 



fields (/> which are periodic in the and satisfy Dirichlet boundary-conditions at = 0, L. 
The first thing to do is to perform the gaussian integrals to get (A' is the Laplacian taking 
derivatives w.r.t. x') 



2_ Of Of 

D^{x^) = exp{—K{x,x)yexp{±--JA'K{x,x')d^x'±--jA'K{x',x)d^x'} (54) 
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where the integration is over x' and the kernel (employing the Schwinger mass (^)) 



is with respect to Dirichlet boundary conditions. From its explicit form one finds |17] 



^(^'^) = ^ E (l - cos(2n7re)) ( ^^^(^^'"^ _ (n ^ ^ n'^ + il^L/T^f )) (56) 

n>l 

/rr. - ^g-O I ch(^L(i-20)-ch(^L) 

^ " 2r + 2/i/5sh(/iL) + 

1 ^ ch(m7r/r) - ch(m7r(l - 2^)/t) . i ; , ^ , „ , 

— > — ^ — — — ^ — 7^ - (m ^ A/m2 + (/i/3/27r)2 ) (57) 

2-K msh(m7r/T ^ V ^'"'"z ; 7 v ; 

m>l ^ ' ' 

which is perfectly finite as well as 

sh(/iL(l-e))+sh(/iLO 



A'i^(x,x') d^x' = / A'i^(x',x) d^x' = ^ — -^^^V^ ^ - 1 (58) 

sh(^L) 



from which the factor (54) can be computed. As one can see, all three expressions (56, 57 



tend to zero as x approaches the boundary. The result for ( |54| ) is then found to read 
D^{x') = exp{^ E (i _cos(2nvrO)(^^^^ - (n - ^n2 + (/.L/vr)2 )) } • 

i n>l 

rl ^ ch(m7r/r) - ch(m7r(l-2^)/r) l~271~7hr\2 \\ 
m sh(mvr/r) (- - V-^ + (M/5/2vr)2 )} • 

exp f ± ^ ."(^M^^(l-e))+sh(M6 _ . I 
I- V sh(;uL) / J 

which is independent of x^ as required by translation-invariance. 

3.3 Analytical Expression for Condensate at Arbitrary Points 



The chiral condensate (|4^) is simply the product of ( |5^ ) and ( p9D or equivalently ( |53D and 
(pO). Rescaled by natural units (H), it takes the form 



(V't P±V')(a;^) g±e.ch(A^(l-20/2)/ch(A^/2) 

(|e|/0F) ^ ^ 4A 

^ sin(7r,^)ch(7rnr) J cos {2TTnc) 9^^ {Cjir) dc 
Sin2(7r0+sh2(7rnr) / <^ (c,ir) dc 

{l^ E (l - cos(2n<)) (^^M!^ _ ^ y'n2+Ar^(A/7r)2)) } (61) 

/ n>l ^ 
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^±0-ch(A/^(l-2O/2)/ch{A^/2) 



2a 



E(-i) 



sh(7r(m+^)/T) 



^ / C(l-g) ch(A77^(l-20)-ch(A7]^) x^ 



(62) 



exp 



r 1 ^ ch(m7r/r) - ch(m7r(l-20/T) , l ' ~ ,^ \o ^^ 



m>l 



where the ^-independent but ^-dependent parts of the two factors have cancehed. Note that 
the two forms ( |6T|) and (|6^) are identical for any finite (dimensionless) lengths a and A, but 
enjoy excellent convergence properties in the low- (r 1) and high- (r <C 1) temperature 
regimes, respectively. 



4 Numerical Evaluation for A^/ = l and Nf = 2 



The general result (61 , |62| ) shall be specialized to the cases Nj = 1 and Nf = 2 as the integration 
over c can be done in these cases. The aim is to collect specific observations concerning the 
difference between the single-flavour and multi-flavour case. 



4.1 Specialization to = 1 at Arbitrary Points 

For one flavour the c-integrals in (^ |6^ ) can be performed. Taking into account the fact 
that they extend over the interval [-1/2,1/2] the result is found to read 



i\e\/V^) 



{\e\/V^) 



„±6»-ch(A/2-A5)/ch(A/2) 

^- E(-i 



sin(7r^)ch(7rnT) 



4A 



exp I ^ 1^1 — cos(2n7r,^)^ ^ 



cth(n7rT^ 



exp(— n^vrr) 



n>l 

„±e-ch(A/2-A5)/ch(A/2) 



n 



±- 



n 



1 



(A/vr)2))} 



2(7 



g7r(m+0(2fc+m+C)/r(-gj.£(- 

fcez 



(fc+m+g+l/2)v/¥ N ^^^. (fc+m+g-l/2)^/¥ ^^ 
\/t ' ^ '''' 



f ^ g(l-0 , ch(A(l-2g))-ch(A) ^^ 

r(— ^ + — ^^hm — ) I • 



exp 



m>l 



h(A) 

ch(m7r/r) — ch(m7r(l — 2,^)/t) 
m sh(m7r/T) 



(63) 



(m ^ y^m2 + (a/27r)2)} (64) 



where again the two equivalent representations (|6^) and (^) enjoy excellent convergence 
properties in the low- (r 1) and high- (r <C 1) temperature regimes, respectively. 
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4.2 Specialization to = 2 at Arbitrary Points 



For two flavours the c-integrals in ( |6l| , |6^ ) can be performed. Taking into account the fact 
that they extend over the interval [-1/2,1/2] the result is found to read 

iij^Piij) _ e^0.cH\/V2-V2XO/cHX/V2) sin«)ch(7rnr) kez 



±- y (-1 



(|e|/V^) 4A 'r'y sin2«) +sh2(^nT) ^ e"2fe2^- 

kez 

exp {i ^ (l - cos(2nvrO) (^^Ml^ _ ^ y'n2 + 2(A/^)2 )) } (65) 

n>l 

{^^P±iP) ^±e-cHx/V2-V2XO/cHx/V2) i 
i\e\/^) - ^ 2a sh(7r(m + O/r) 

^ ^ l+(-ir+'^ ^,r((p+m+f)^-P^-q^)/2r(-^j.£(- (p+m+g+l)v^ -) ^^^^^ (p+n^+g-l)0F ^^ 



2 E e-'^92/2T 
ex r^/ e(l-0 ch(V2A(l-20)-ch(V2A) x^ 



2 V r V2cT sh(^/2A) 

fl ^-^ chfmvr/r) — ch(m7r(l — 20/t) , I ' ; TT, ~i , , 

l2 mshfmvr/r ^ V w y ;j \ i 



m>l 



where again the two equivalent representations (65) and (pq) enjoy excellent convergence 



properties in the low- (r ^ 1) and high- (r ^ 1) temperature regimes, respectively. 
4.3 Numerical Evaluation of (^-Dependence 



We are now in a position to evaluate formulas (^^) for Nf = l as well as formulas (55,pq) 
for Nf = 2. The first hint regarding the difference between single-flavour and multi-flavour 
cases might come from observing how the spatial dependence of the condensate behaves as 
temperature and box- length vary. Since the condensate diverges on the boundaries as 
and (1 - ^)-^ respectively Q, it is the quantity 4^(1 - ^)| (^/;tP±^/;) |/(|e|/y?) at 6^ = which is 
displayed in figures I and |. 

From inspecting figures ^ and |^ one is lead to the observation that for any finite tempera- 
ture the condensate ends up "creeping into the boundaries" (and fading away at any internal 
point of the box) once the box-length is sufficiently large. 

On the other hand, cooling the system seems to have the opposite effect, i.e. to a first 
approximation, the shape of the spatial distribution of the condensate seems to depend on the 
ratio T = a/2X. However, as a and A both tend to be large, the spatial distribution starts to 
look qualitatively different in the two cases Nf = l and Nf = 2, respectively. This is a pictorial 
hint indicating a difference between the single-flavour and the multi-flavour versions of the 
model once the twofold limit a — > oo, A — > oo is performed. 

^ There is no reason to start worrying: On the boundary the field has to satisfy a boundary condition which 
lets the free Green's function take the form given in ([l9| ) . This expression stays finite as one of its two entries 
- X OT y - reaches the boundary but not if x and y reach the same point on the boundary, which simply means 
that there is no consistent double-boundary solution. 
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4.4 Specialization to Nf = 1 at Midpoints 

As one concentrates on the midpoints = 1/2), formulae (]63|) and take the form 



(|e|/0F) 47r V ' ch(n7rT) 



exp{7-2^(-iyi^o(jA)} 



n>0 



(|e|/0F) 47r ' sh(7r(2m + l)/2T) 

^ch( -(^"-+y+^^ )(erf((^±^)-erf((^)) 



fe>0 

7r(l -th(A/2)) 



^ i>i 

-P { - 2 E J/. , n - (-^v'-^+W2vr)2)} (68) 

m>l ^ ' 



respectively (for details, see appendix), from which we derive 

{^'^P+^)(-\ g±e/ch(A/2) f ^ . . 

lim ^ , ifi^,^^ = exp{7-2^(-l)^KoO-A)} (69) 



hm = ±- 

A-.c« (|e|/^) 47r ^.^^ 



i-exp{7-2^i^o(j<T)} (70) 

for Nf = l. Thus 



A^oocr-^oo (|e|/A/vr) 47r a-^oo x^oo {\e\/^/^T) 

for Nf = l. Note that ( |70D and the second equality in ( [tTD correct for an erroneous result in 
1 17 1 which was won by performing the limit under the c-integral. 



4.5 Specialization to A(/^ = 2 at Midpoints 

As one concentrates on the midpoints = 1/2), formulae ( |65| ) and (|6^) take the form 



(|e|/v^) 40rvr 

-n'^TTT/2 



ch(n7rr) ' 1 + 2^6- 



i + 2E(-i 



fc>l 
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n>0 



" ^ 4VFV^ 4^/"^^ sHn{2m + l)/2r) ^ 

qeZ p>0 ^ V2^/'^ V2^/'^ 

q£Z p>0 ^ V2r/7r VW'r 



r7 7r(l-th(A/V2)) / 



respectively (for details, see appendix), from which we derive 



^^^^^Wr - " 4Vvr exp{|-E(-l)^i.o(.V2A)} (74) 



hm ^ = (75) 

A-.^ (|e|/V^) 

for iV/ = 2. Thus 

lim lim 7*^/^)1^ = = lim lim r'^^)!^'^ (76) 
A^ooo-^oo (|e|/v7r) (T-+ooA^oo (|e|/V7r) 

for Nf = 2. Note that in the first case the condensate decays rather reluctantly (oc 1/\/A) 
only under the outer limit, whereas in the second case it decays exponentially fast under the 
inner limit already (see appendix). This difference may be stated in equations through 

li„ ,i,„ = (78) 

o-^ooA-*oo (|e|/0r) 

which is peculiar for Nf = 2. 

4.6 Numerical Evaluation at Midpoints 

We are now in a position to numerically evaluate formulas (67, ^) for = 1 as well as 
formulas (|7|, |7|) for Nf = 2. 
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Figure 4: Hip'' P±ip)\/ m as a function of log{kT/fii) at L = O.OS/m for Nf = 1 and iV/ = 2. 
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Figure 5: Hi^^ P±ip)\/ m as a function of log{kT/fii) at L = 27.0/^i for Nf = 1 and Nf = 2. 



The following figures display the absolute value of the condensate at = in the center 
of the box. The two-dimensional graphs show the condensate as a function of (T = /3 • |e|/7r^/^ 
(or a function thereof) at fixed value of A = L • |e|/7r^/^. Throughout we use /Ui = |e|/7r-^/^. 
The surface- and density-plots show the condensate versus (a function of) a and A. 

It might be worth mentioning that in each of those figures both representations - ( ^2[ ) 
and ( |84|) in case of = 1 as well (^) and ( |87|) in case of = 2 - were used, the former 
formulas got evaluated for low temperatures (r 1), the latter ones got evaluated for high 
temperatures (r ^ 1) — the switching being done in the region of the crossover transition, 
which essentially gives a numerical check that the representations (^) and (|8^ as well as the 
representations (36|) and ( pTj ) indeed might be identical. 

and |5| the system is seen to undergo a surprisingly well-localized crossover, 
- doesn't meet the criteria for a phase-transition of whatever kind as it is 
arbitrarily smooth (i.e. G C°°). 

By comparing figures ^ and |5|, one realizes that increasing the box-length essentially moves 
the "kink" to the left, i.e. increasing the box-length results in a decrease of the "critical 
temperature" (the effect being much stronger in case Nf = 2 than for Nf = 1) and one starts 
wondering whether this kink-phenomenon survives the limit L — > oo. We will see that the 
answer to this question depends in a critical way on the number of flavours; the two cases 
Nf = 1 and Nf = 2 turn out to be different. The "plateaus" seem to be equally high for one 
and two flavours as long as the box-length is small (cf. figure §) but the height of the plateau 
decreases unequally rapidly if the box-length is increased (cf. figure - compare the scales !). 



From figures 
which - however 
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Figure 6: |(V'^-f± "0)1/^1 ^ function of log{kT / fii) and log(L/ii) for Nf = 1 and Nj = 2. 

In order to reach an intuitive understanding it is worth having a look at figure |^. Thereby 
one reahzes that there is a considerable fraction of the log(A;T//ii)-log(L/xi)-plane where the 
condensate is exponentially close (but not equal) to zero. Thus it might be helpful to introduce 
the concept of a quasi-phase with almost restored chiral symmetry. In order to do so one 
has to decide on a trigger-value which the condensate has to exceed in order to constitute 
a point (in the log{kT/^i)- log(L/ii)-plane) with manifestly broken symmetry. While this 
choice is - in principle - arbitrary, it seems natural to agree on half of the classical value 
of the condensate in the original (one-flavour) Schwinger model as the discriminator which 
makes the distinction between the two "quasi-phases" . Numerically, it is about 0.07. 

Doing so results in generating the two contour-plots in figure |^: White points are those 
which satisfy the criterion \ {ip'^P±ip)\ > e'^^/e-^/vr/Svr - they constitute the (quasi-)phase with 
manifestly broken symmetry. Black points are those which satisfy the criterion |(V'^-P±'0)| 
< e"^ ^/e^Jn/Sn ~ they constitute the quasi-phase with quasi-restored symmetry. These two 
quasi-phases are separated by a "crossover-line" which essentially corresponds to the "border" 
of the zero-level plane in figure In the form shown in figure ^ the concept of quasi-phases 
proves useful as it clearly shows in which areas of parameterspace the cases Nf = l and Nf = 2 
seem to be similar and in which areas each of them shows a clearly distinct behaviour. 

First we notice that the point (3 = L = (lower right corner in figure ^) seems to lie close 
to or right on the crossover-line in either case = 1 and Nj = 2 (note that data near the 
boundaries are cut off for numerical reasons). This observation lets us go back to formulas 
(|6l| ) and (62) and derive the (moderately interesting) noncommutativity-phenomenon 



Jhn lim (V'tp±V')(^) = (ViV^) (79) 

Ihnlhn (V^tp^V'X^) = oo (ViV;) (80) 

which is universal for any number of flavours. 

Second we notice that the point /3 = L = oo (upper left corner in figure |^) definitely belongs 
to the manifestly broken (pseudo-) phase for Nf = l (l.h.s. of figure but for Nf = 2 (r.h.s. of 
figure]^ the point f3 = L = oo seems either to be part of the quasi-phase with almost restored 
symmetry or to lie right on the crossover-line (note, again, that areas close to the boundaries 
in figure ^ are cut off for numerical reasons). This observation lets us try to zoom into the 
upper left corner in these quasi-phase-structure plots - an attempt which results in figure @. 
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Figure 7: Quasi-Phasestructure as a function of log{kT/fj,i) and log(L/ii) for = 1 and 
Nf = 2, respectively. Definition: see text. 



From the l.h.s. of figure ^ we learn that the one-flavour-system approaches its standard- 
value for the order-parameter in a very unspectacular way: The "classical" value for the 
condensate in the single-flavour model represents a plateau which is reached smoothly from 
any side. From the r.h.s. of figure |8| we learn that the two-flavour-system behaves in the large- 
L-limit in a way which depends rather sensitively on whether the temperature is exactly zero 
or finite: For any finite temperature the limiting value (zero) is approached smoothly (a state- 
ment which extends to the first derivative of the condensate), whereas at zero temperature 
the condensate seems to display a "square-root type" behaviour (as a function of 1/A). This 
is exactly what was predicted by formulas (^ 75): At zero temperature the two-flavour 



condensate goes to zero under A ^ oo - but only very reluctantly: oc 1/\/A (up to exponen- 
tially small corrections). On the other hand, for any fixed finite temperature, both the the 
condensate and its first derivative w.r.t. 1/A vanish exponentially fast under A oo. 



4.7 The order of the phase-transition at Tc = 

The statement that in the two-flavour model the behaviour of the chiral condensate as a 
function of A = /iiL depends in a very sensitive way on whether cr = /ii/3 is finite or infinite 
turns out to be so crucial in the following that one would particularly welcome some further 
numerical evidence that there is, indeed, a "square-root type" behaviour of the condensate 
for 1/A^l at 1/(T = (as opposed to a smooth behaviour for A — > oo at l/cr>0). 

Though there is, quite generally, no numerical proof for smoothness, strong numerical evi- 
dence can be given that the boundary associated to T = of the two-flavour condensate surface 
in the r.h.s. of figure ^ does not just look like a square-root but, indeed, asymptotically gets 
a square-root and that this small- 1/A-behaviour is, indeed, specific for r = 0. To this end we 
simply decide to plot, for the two-flavour system, the quantity I'^Tx^l'^e'^ -^fX |(V'^-P±V')l//^i) 
i.e. to include a factor cx The result is shown in the r.h.s. of figure]^: the r = 0-boundary 
tends to 1 whereas the 1/L = oo-boundary seems to be compatible with 0. In summary, the 
r.h.s. of figure ^ provides an independent numerical check of the analytical work (presented 
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Figure 8: Zoom-out of the area C7 3> 1, A» 1 (upper left corner in ^ for = 1 and Nf = 2. 

in the appendix) which has been done to get from ( [72| , [7^ ) to ( [7^ ^). 

At this point a pitfall arizes: If, for the two-flavour case, the "alternative quantity" 
2V47]-i/2g-7 . P±ip)\/ f^i would be an order-parameter (actually it is not an order- 

parameter, as we shall see in a moment), the very fact that it acquires a finite value (1 in our 
normalization) at P = L = oo if the zero-temperature- limit is taken prior to the limit of infinite 
box-length while staying zero if the two limits are performed in reverse order would lead us 
to the (false) conclusion that at P = L = oo the system consists of two coexisting phases with 
relative weights which could be set by the way the point P = L = oo is approached: Following 
the corresponding "line of constant altitude" in the r.h.s. of figure the relative weight of 
either phase (seemingly) could be given an arbitrary value between and 1. This means 
that we would reach the (incorrect) conclusion that the two-flavour Schwinger model shows 
a first-order phase-transition with critical temperature Tc = 0. 

While the statement that the two-flavour Schwinger model exhibits a phase-transition at 
Tc = is indeed correct (see below), a classiflcation as of first order is not - the transition is 
actually second order. To this end we shall first point out the fallacy in the above pseudo- 
reasoning and then present the correct argument. 

The first question is: what is wrong in the argument which was declared to be a pitfall ? 
The answer is: nothing is wrong - except for the fact that the if-condition is not fulfilled. 
The point is: the "alternative quantity" (the condensate times a factor \/A) is not an order- 
parameter. This might come as a surprise: Usually, in order to study critical phenomena or 
spontaneous symmetry breaking in field theory, the prescription is that one should choose 
an operator which is non-invariant under the symmetry in question. Then, in the presence 
of an external symmetry-breaking field, the vacuum-expectation-value of this operator for 
different temperatures and different external fields indicates the type of the phase-transition. 
In the case at hand the true order-parameter (the chiral condensate) and the "alternative 
quantity" differ from each other by a scalar factor, i.e. these two operators transform in 
the same way under a chiral transformation. Nevertheless, the correct order-parameter and 
the "alternative quantity" do not prove equally useful to study the chiral structure of the 
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Figure 9: Left: The order-parameter KV'^-Pi^)!//^! of the one-flavour system as a function 
of the dimensionless temperature at infinite box-length. Right: Zoom-out of the upper left 
corner in the r.h.s. of figure for the two-flavour system where the alternative quantity {not 
an order-parameter !) 2'^^^TT^^'^e~'^ ■ ^^^'^\{ip^ P±ip)\/ Hi is plotted rather than {ip^ P±ip) / Hi - 
providing direct numerical evidence for the noncommutativity phenomenon [7^ ). 

theory: The "alternative quantity" can be seen as a specific member of a one-dimensional 
class of operators where the parameter is just the exponent of the prefactor A in front of the 
condensate. Operators for which this parameter is smaller than 1/2 tend to zero under A — > oo 
whereas operators for which this parameter is bigger than 1/2 diverge under A ^ oo (at zero 
temperature). Demanding this parameter to equal 1/2 produces the finite jump at critical 
temperature shown in figure |^ but there is no physical meaning in this behaviour whatsoever: 
The "alternative quantity" oc A^/^- P±ip)\/ Hi can not be used as an order parameter since 
its very definition either requires analytical knowledge about the behaviour of the condensate 
(which is the true order parameter) for large values of A at zero temperature or so-to-say 
"critical tuning" which is completely intolerable. Moreover, even if we had overlooked this, the 
(incorrect) conclusion that the phase-transition is first order would not even be self-consistent: 
Whenever a system exhibits a true first-order transition a sufficiently small perturbation by 
the symmetry-breaking external field still leaves the transition first-order. This, however, is 
definitely not true for the case at hand: The condensate (and henceforth any operator in the 
one-dimensional set introduced above) was found to show an arbitrarily smooth crossover (as 
a function of a) for any nonzero value of 1/A (which, in our setting, plays the role of a small 
quark- mass term). 

Having convinced ourselves that the condensate P±ip)\/ Hi is the unique legitimate 
order parameter it is straightforward to read off the order of the phase-transition of the two- 
flavour system directly from the r.h.s. of figure g At zero temperature the condensate follows 
a "square-root" behaviour down to its limiting value (sc. 0) under 1/A ^ 0-|- whereas for any 
finite temperature it decays exponentially and henceforth such that even its derivative tends 
to zero under 1/A — > 0-|-. Thus there is no finite jump and no phase-coexistence at Tc = 0, the 
susceptibility (which, in our setting, is associated with the derivative of the condensate w.r.t. 
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1/A) diverges under 1/A ^ 0+ for r = like 1/V(1/A) but tends to zero under 1/A — > 0+ 
for T>0 and the transition is second order. 

4.8 Determination of d 

A system with a second order phase-transition exhibits a critical behaviour at its vicinity 
which is described by a set of critical exponents. While some of them refer to the broken 
phase and thus can not be defined in the case at hand, some other describe how the system 
behaves at or slightly above the transition under the influence of a symmetry-breaking external 
field and may be defined even if Tc = 0. 

In order to rephrase our results in the language of the usual approach where the infrared 
regulating and explicit chiral symmetry breaking device is the small fermion mass m rather 
than the inverse box-length 1/L one has to agree on how one of these two quantities shall 
be translated into the other one. It is clear that 1/L must be associated with a monotonic 
function of m. Sticking to the "naive" (i.e. dimensionally motivated) choice of identifying 
1/L <-> m, the order-parameter in the massive zero-temperature Schwinger model is easily 
derived from ( p9| ) (for Nj > 2) 

hrn (V'^PiV) oc iv//'"^^ mW-i)/^/ exp {^(7 - 2 ^(-ly i^o(j V^^))} (81) 

whereas at finite temperature the order parameter and its first derivative w.r.t. m still go to 
zero under m 0+. From this we learn two things: 

(1) The phenomenon of a second-order phase-transition with zero critical temperature is 
not specific for the two-flavour model - the (massless) Schwinger model has a phase-transition 
at Tc = for any number of flavours bigger than one. 

(2) The critical exponent 6 which is defined through {ijj'^P±ijj) oc m}/^ (at T = Tc) is 
Nf /{Nf—1), i.e. 6 = 2 for Nf = 2 (upon identifying 1/L ^ m - for a critical remark see below). 

5 Discussion and Conclusion 

The present paper has been devoted to a study of the A^-fiavour euclidean Schwinger model 
on a finite-temperature cylinder with SU{Nf)A breaking local boundary conditions at the two 
spatial ends. We have investigated the value of the dimensionless condensate {tp^ P±ip) / fJ-i at 
midpoints ~ which we used as order-parameter - as a function of the dimensionless inverse 
temperature a = iii(3 and box-length X = fiiL. 

Our aim was to give a qualitative picture of the behaviour of the Schwinger model when 
quantized as described above. We found that on a logarithmic temperature-scale (with L kept 
finite but fixed) the condensate undergoes a well-localized crossover from a fairly constant 
value (at low temperatures) to a value which is exponentially close to zero at sufficiently high 
temperature - both for Nf = l and Nf = 2. From this we concluded that it is most reasonable 
to distinguish a (quasi-) phase with manifestly broken chiral symmetry from a quasi-phase 
where the chiral symmetry is almost restored - the distinction being done trough a critical 
value of the condensate which was defined as half of the classical value of the condensate in 
the single-flavour model at zero temperature. 

The numerical illustrations of our analytical results show that there is a qualitative dif- 
ference in the behaviour of the single-flavour model as compared to the two-flavour model if 
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the combined limit of large box-length and small temperature is considered. For Nf = l the 
condensate takes a nonzero (T-dependent) value after L ^ oo. For Nf = 2 the condensate 
ends up being exactly zero under L ^ oo, but the way it approaches this limiting value is 
rather different for fixed T = versus fixed T > 0: extremely reluctantly (i.e. oc 1/^/L) in 
the first case versus exponentially fast in the second case. This means that the susceptibility 
(defined as the first derivative of the condensate w.r.t. 1/L) diverges in the two-flavour model 
like \/L under L ^ oo at zero temperature but tends to zero under L ^ oo at any fixed 
finite temperature. This means that the two-flavour system exhibits a second-order phase- 
transition with critical temperature Tc = - a result which was shown to extend to the cases 
Nf > 3. For = 1, on the other hand, the order-parameter stays nonzero after L ^ oo 
for any T, and the fact that it is just exponentially close but not equal to zero for T ^ |e| 
illustrates the statement by Dolan and Jackiw that the one-flavour system never restores the 
anomalously broken C/(l)^-symmetry at finite temperature |^. 

It is worth emphasizing that the results summarized so far provide direct pictorial evidence 
for the claim by Smilga and Verbaarschot which was based on an interesting indirect 
argument: These two authors used the result for the scalar susceptibility in the Schwinger 
model with degenerate massive flavours determined via bosonization rules and completed the 
list of critical exponents at small but nonzero temperature. Their result was that this list 
can be understood most easily as to consist of entries which satisfy the scaling relations for a 
system slightly above a second-order phase-transition with Tc = 0. 

There is, however, one important point of numerical disagreement: Based on our analytical 
result plus the identification rule m <-> 1/L we have found the critical exponent 5 = 2 for the 
two-flavour system - which is at variance with the bosonization-rule based result 6 = 3 (which 
also agrees with the mean- field value) found in the literature 

As far as our part is concerned all we can say is that we have tried most diligently to 
make sure that our result is not due to an error in the analytical computations: Our finding 
for the numerical value of 6 in the two- flavour case stems from (74) which was derived from 
(|7^ ) . Through an explicit plot (figure |9|) we tried to convince ourselves numerically that all 
the subleading factors in ( [7^ ) do indeed get marginal in the zero-temperature limit and the 
overall A-dependence tends to (for Nf = 2). Moreover numerical evaluations of ( |72|) and 

( [73| ) (after having cut down the infinite sums to an appropriate finite number of terms) were 
found to agree to more than a dozen decimal places for a large variety of a- and A- values. 

As far as the result for 5 in the literature is concerned all we can say is that the results 
for the massive Schwinger model we are aware of are based either on perturbation-theory in 
the fermion-mass or on bosonization-rules (see |25] and references therein). While the 
original derivation of some of these rules was again within the mass-perturbation approach 
[p6|| , the rules themselves seem not to be tied to perturbation-theory: Hetrick, Hosotani 
and Iso found in the bosonized massive Schwinger model with 2 or 3 degenerate flavours 
a noncommutativity-phenomenon between m — > and L — > oo (where L plays in their 
scheme the role of a finite inverse temperature (3) [p7|], which would be rather surprising if 
the bosonization-rules would not go beyond mass-perturbation theory. Nevertheless, for the 
bosonization rules to be applicable the condition m <C |e| has to hold true, whereas our results 
stem from an approximation-free analytical computation and are supposed to be exact for 
any inverse box-length and temperature. 

At this point it should be stressed that the identification 1/A <-> m (which was necessary 
in order to extract, from our results, a value for 5 in the usual approach where the symmetry 
breaking field is the quark-mass rather then the inverse box-length) is not canonical; a finite- 
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volume effect might, in principle, spoil the validity of our simple dimesionally motivated 
identification rule. Though this possibility seems rather unlikely, we feel that in the light 
of the known peculiarities of the two-dimensional world (the classical one being described 
in 0) it can't be ruled out a priori. Needless to say that in the present situation both 
further analytical results in the massive multi-flavour Schwinger model and a determination 
of its critical indices from the lattice would be highly desirable. However, the fact that the 
critical temperature is zero provides a sort of a challenge for the lattice approach as it requires 
sophisticated reweighting techniques. 

In summary we have presented a study of the chiral condensate in the one- and two- 
flavour finite-temperature Schwinger model in a quantization scheme where the usual quark- 
mass term is replaced by bag-inspired chiral symmetry breaking boundary conditions. Un- 
like results won from bosonization-rules or within mass-perturbation our formulas represent 
(hopefully) exact analytical findings gained through a straightforward evaluation of the path- 
integral. We have introduced the concept of quasi-phases in order to distinguish regions in 
the parameter-space where the symmetry under investigation is manifestly broken from those 
where it is almost restored and we have provided a direct pictorial verification of the claim by 
Smilga and Verbaarschot that the two-flavour Schwinger model undergoes a phase-transition 
at Tc = and that the transition is of second order . At the time being we are unable to 
resolve the discrepancy between our value for the critical exponent 6 = Nf/{Nf — l) and the 
bosonization-rule based result 6 = {Nf + l)/{Nf — l). 
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Appendix 



Here we shall give the details for the derivation of formulas 



As one concentrates on the midpoints (^ = 1/2), the r.h.s. of (63) takes the form 

g±e/ch(A/2) 



(63) 



4A 



exp 



cth((2n-M)7rT) 

±9/ch(A/2) 



n>0 



(2n + l) 



((2n+l)^^(2n+l)2 + (A/vr)2)} 



4A 



(l + 2i:(-ir5;(^exp(-„V)). 

1 



n>l 



exp 



I (2n+l) 



„>o(2n+l) V(2n+l)^ + (A/vr)^ 

g-2'm(2nfl)7rT 



■■} 



exp 



{4 E E ^2n + l) - ((2^ + l)-V(2^ + l)' + (Vvr)2)} 



n>0 m> 
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±e/ch(A/2) . ^ 1 

±- (l + 2y(-l)" — expC-nVr 

n>l ^ ' 

n 1, A 
r 1 

exp 



exp {2(1 + ^ log(^) - Y.{-iyK,{3\)) ] 

£ (2„ + l)(e«U.^-l) - (P»+l)-v'P» + l)' + (W)} 



n>0 

which is the result (|67|) quoted in subsection 4.5 and where we used |22] 



with 7 = 0.57 . . . the Euler gamma and Kq the modified Bessel function. 

As one concentrates on the midpoints (.^ = 1/2), the r.h.s. of (64) takes the form 



g±0/ch(A/2) 



^^^^ ^ 2CT «h(vr(2m + l)/2r) 



meZ 

^ g7r(m+l/2)(2fc+m+l/2)/T (^gj.f (- {^±!Zl±l)v^^ _gpf |- (fc±!B)v^^-) 

fcez ^ 

2 

f / 1 th(A/2)N ^ th(m7r/2r) / n,^ /o ^2^\ 

expj^i^ + (m^^m2+(a/2^)2)| 

ra>l 

±0/ch(A/2) ;l 



2a ^^^^ sh(7r(2m + l)/2r) 



^ g7r(m+l/2)(2fc+l-m-l/2)/T(^gj.£(^{l+(^+l))vArs| | erf ( ^^"^^^^ ^^^ )) 



-2n(m7r/2T) ^ 

e^p{2E E(-l)" (m-y'm2+(a/27r)2)} 

m>l n>l 

„±e/ch{A/2) „-7r((2m+l)2-l)/4r 

a ' sh(^(2m + l)/2r) 



fc>0 



m>0 

^ eh( -(^-+^if^-+i) )(erf( »^^+^^>tW- 



r , vr 1-th A/2 

exp 7 + log — + ^ - 2 ^ Ko jcT 

47r a 

{ - ^ E J/. ^ 1) - (n.^ V^^+W2vr)^ )} 



m>l 
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which is the result (|6^) quoted in subsection 4.5 and where we used |22] 

with 7 = 0.57 . . . the Euler gamma and Kq the modified Bessel function. 

As one concentrates on the midpoints = 1/2), the r.h.s. of (65) takes the form 

E„ — fc^TTT— (n— fc)^7rT 
^ ^ ^^^^ 

= ^ iA ch(n7rT) E 

cth((2n + l)7rr) 

exp ' 



-r?-KT/2 g-(n/2-fc)-'27rT_)_g-(n/2+fc)-' 27rr 

y(-i)" 



4A ^ ch(n7rT) " ^ 6-2^^ 



''''P{^j2n+1) " V(2n+l)2+2(A/7r)2 } 

-2m(2r)+l)7rr , 

exp{2 5: y^^-^-((2n+l)-y(2n+l)2 + 2(A/vr)2 



n>Om>l 



^l + 2y(-l)" 



4A V ^ ch(n7rr) 1 + 2 e-^^'^'^^ 

~ fc>i 



2 + 2^°^^^' 

{'^. (2n+l)(e2(U^--l) -(^'"+^^^V('"+^)' + '(^/")'^} 



exp {| + - log(^) - Y.i-^yK,{j^ A)} 
r 1 

exp 



n>0 



which is the result (^) quoted in subsection 4.5. and where we used (pS]). 

As one concentrates on the midpoints = 1/2), the r.h.s. of (66) takes the form 

„±6l/ch(A/V2) 1 

(66) = ±- y (-1)™ — — — X 

^ ' 2cT > sh(7r(2m+l)/2T) 

^ ^ ^^(fo+^+i/2)2-p2-q2i/2r^^^f^ (p+m+3/2) ^ ^^^^ (p+n^-1/2) 

2 I] Q-T^q'^/'iT 

rvr/l th(A/\/2)\ 1 ^ th(m7r/2r) , / o /o n9n1 
exp { - (- - + - E - (- - yrn2 + 2(./2vr)2 )} 

^ m>l 
„±6»/ch(A/v^) p-7r(m+l/2)2/2T 

2a ^ sh(vr(2m + l)/2T) 
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J2 g-7rg2/2r^ c^fo+l/2)(m+l/2)/r (^^f ( l+(p+l/2) ^ erf ( ^^^^^^^^^ )) 

qeZ peZ ^ V^r/TT y/2T/TT 

fTT/l th{X/V2) \ 1 V- J_ 1 -1 

""""PUUr" V2a ^ + 2 m " + 2(,j/2^)2 / 

-2n{m7r/2T) ^ 

exp{ E E(-ir {m^^m^ + 2{a/27ry )} 

m>l n>l 

„±e/ch(A/v^) „-7r((2m+l)2-l)/8T 

2a sh(7r(2m + l)/2T) 

J2 e--3V2r ^ ^^^( x(p+l/2)(m+l/2) ( (p+l/2) + l ) ( (p+l/2)-l 

^ g-7rgV2r ^ ( l)''+''-"'sh( ) (erf ( + ^ ) ci'f ( )) 

geZ p>0 ^ \/2^ 

2 J2 e-'^9V2-r 

7rth(A/\/2) 1/ vr ^ ,\/2"cj, 

i>i 

r ^ 1 . / . ^, .o . 

exp 



f 7rth(A/\/2) 1/ vr , ,V2cr, , ^ AT 

{ - E ^ , r '"^ v^m^ + 2(^/2,)'2 )} (87) 



(gTrm/r _^ 

which is the result (^) quoted in subsection 4.5. and where we used ( p5| ) 
For later use we finally mention that 



-7r((2m+l)2-l)/8r 

X 



sh(7r(2m + l)/2r) 
j2 e-^gV2r ^ ch( -(^+^/^»"'+^/^) )(erf ( ) -erf ( ))+ 

J2 e-'f'?V2r J2 ( l)P+q-mg|^/ ^(p+l/2)(m+l/2) (p+l/2) + l x ^^^z (p+l/2)-l 

ggz p^o ^ 

~ 4 * exp(— — ) 
4r 

for T <SC\ and that the limit r — > oo for arbitrary takes the form 
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